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Abstract 

We employ a simple effective model to study the chiral dynamics 
of two flavors of quarks at finite temperature and density. In particu- 
lar, we determine the phase diagram in the plane of temperature and 
baryon chemical potential as a function of the pion mass. An inter- 
esting phase structure occurs which results in zero, one or two critical 
points depending on the value of the vacuum pion mass. 
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1 Introduction 



The physical pion mass is small but not zero. In consequence, the conven- 
tional wisdom is that there is no true thermodynamic chiral phase transition 
at finite temperature T and zero baryon chemical potential fi, at least when 
the effects of the strange quark are neglected. Instead, there is expected to 
be a curve of first order phase transition in the n — T plane which termi- 
nates in a second order phase transition at some critical point {nc,Tc). The 
location of the critical point obviously depends on the physical (vacuum) 
pion mass. This topic has been under intense theoretical study using various 
effective field theory models, such as the Namu Jona-Lasinio model [H [21 E], 
a composite operator model [1], a random matrix model [5], a linear sigma 
model [3], an effective potential model [6], and a hadronic bootstrap model 
[7], as well as various implementations of lattice QCD [HI [9l [TOl [11]. Reviews 
of the subject were presented by Stephanov in the last few years [I2]. It 
is also of great interest because collisions between heavy nuclei at medium 
to high energy, such as at FAIR (Facility for Antiproton and Ion Research), 
which is under construction, may provide experimental information on the 
phase diagram in the vicinity of a critical point. In contrast, current experi- 
ments at RHIC (Relativistic Heavy Ion Collider) and future experiments at 
LHC (Large Hadron Collider) may create too much entropy and therefore 
miss the high baryon densities needed to explore that region. However, that 
depends on the location of the critical point. In this paper we study the 
phase diagram of the linear sigma model coupled to two flavors of identi- 
cal mass quarks. Whereas the sigma model is an oft-used effective model 
that represents some of the essential features of the chiral dynamics of QCD, 
the reason to couple the fields to quarks is less obvious. Why not couple 
to nucleons, delta resonances, etc. instead? One argument is based on the 
existence of the critical point itself. If a critical point exists, then one can 
go around it without crossing the curve of first order phase transition. The 
effective degrees of freedom should not change too much in following such a 
path. Therefore, if quarks are considered to be reasonably useful degrees of 
freedom on the higher temperature side then they should be useful on the 
lower temperature side too. One could make the same argument for using 
nucleons, delta resonances, and so on. Presumably one could get the same 
answers. However, it would entail including all the baryons in the Particle 
Data Tables and a multitude of coupling constants. For this reason we use 
quarks to carry the baryon number but acknowledge the resulting uncertainty 
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in the results. 

The present work extends that of reference [3] in two significant ways. 
First, reference [3] studies the hnear sigma model with quarks but only in the 
mean field approximation. The present work goes beyond that by including 
thermal fluctuations of the meson and fermion fields, which can be impor- 
tant at finite temperature when the magnitude of the fluctuations becomes 
comparable to or greater than the mean values. The technique we use was 
developed in several papers [TSj [m [15] . Second, reference [3] fixed the pion 
mass at its vacuum physical value, whereas the present work does a scan of 
pion masses from zero to over 300 MeV, which may be particularly useful 
for comparison with lattice gauge theory with different quark masses. The 
present work also extends that of reference [TB] to nonzero chemical poten- 
tial, otherwise the techniques are the same. Even when all other parameters 
of the model are fixed, we find an interesting phase diagram which may have 
zero, one or two critical points depending on the value of the vacuum pion 
mass. 



2 The Linear Sigma Model with Quarks 

The Lagrangian is 

l^ = \ {d^-^f + \ {d^^f - U[cy, tt) + ^ - (a + %^^r ■ tt)] ^ (1) 

where 

[/(a,7r) = ^(a^ + 7r^-/f -i/a (2) 

in an obvious notation. The SU{2)l x SU {2)r chiral symmetry is explicitly 
broken by the term Ha which gives the pion a mass. The scalar field has 
a nonvanishing vacuum expectation value v determined at the classical level 
by the equation 

\v{v'-f)=H. (3) 

The scalar field is thus represented by a condensate plus a fluctuation a = 
V + A. The quarks have no intrinsic mass and only acquire one because of 
the condensate: niq = gv. The four parameters in the Lagrangian, /, H, A 
and g are constrained by fixing the pion decay constant /,r = 92.4 MeV and 
the a mass = 700 MeV, while the quark mass is set to one-third of the 
nucleon mass rriq = 313 MeV. The last required piece of information is the 



3 



vacuum pion mass m^, which is varied from to m„l2 = 350 MeV. At the 
classical level there are the simple relations 

H = Uml 

mz — ml 



2 ^ ml - 3ml 2 

, = ^. (4) 

All of this is standard practice, with the specific choices of /,r, and mq 
corresponding to [I6]. The value of the pion mass in the latter paper was 
fixed at its observed vacuum value of about 138 MeV. 



3 Self-Consistent Method for Determining the 
Equation of State 

The grand canonical partition function is 

Z = I [da] [dTz] [dip] [dip] exp { f dr f d^x{C + |/iV^7V) \ (5) 



where /? = 1/T is the inverse temperature, /i/3 is the quark chemical poten- 
tial, V is the volume of the system, and r denotes imaginary time. In order 
to create an effective mesonic model, we integrate out the quark degrees of 
freedom in the usual way [H], such that 

In Zquark = In dct D (6) 

where D is the inverse quark propagator. This, of course, still depends on 
the mesonic fields via an effective quark mass m, where 

m' = g\a' + 7.') (7) 

which enters into the quasi-particle quark energy E = y/p'^ + m'^. This allows 
us to rewrite the grand canonical partition function solely in terms of mesonic 
degrees of freedom. It amounts to the effective Lagrangian 

C = ^id,n)' + ^id,af~U,sia,n) (8) 
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where 

T 

Ues{(T, tt) = U (a, tt) - — In Zq^^,i,{a, tt) (9) 

is the effective potential. The latter depends on both T and fj, as well as the 
mesonic fields. 

The four equations of motion that follow from the effective mesonic La- 
grangian eq. ([H]) are 

f)TI 

d'd,a+^ = (10) 
d''d,n, + ^ = (11) 

with i = 1,2, 3. Now we decompose the mesonic fields into a condensate 
and fluctuations, as mentioned earlier. We then write out the equations of 
motion in their formal series expansion in the fields and linearize them by 
making the replacement A" n(A"~^)A where the angular brackets denote 
the ensemble average, see below. We make use of the fact that, in this 
approximation, the sigma and pion field fluctuations are independent of each 
other, that is, (A^vr^) = (A")(7r[). The condition for the condensate is 

The quasi-particle dispersion relations for the mesons, namely E"^ = p"^ + 
and E"^ = p"^ + , involve effective masses given by 

"^- = \W/' ^^^^ 

The thermodynamic potential is then computed from 

^ = (f/cff) - lml{A^) - ^ml{7v^) + (15) 

where Q^j and fi^ are the independent particle contributions from the sigma 
and pion quasi-particles. Explicitly 



T 



27r2 



I dpp" -/3E^ + ln(l - e^'^^'^) 



(16) 
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with a similar expression for the pions. 

In this paper we will drop the shift in the zero point energy for simplicity. 
Phenomenologically it may not make sense to keep it since it corresponds to 
high energy/short distance contributions which are probably not described 
well in this effective model. The nontrivial issue of regularization and renor- 
malization in self-consistent approximations has been discussed for the linear 
sigma model several times [THl [13 ED] • 

Ensemble averaging is performed via the technique used in [151 [IB] , which 
in turn is based on the analysis in [T3]; see also [[13]. Consider some arbitrary 
functional tt^), such as lnZquark(o', tt). Expand this functional about 
the average values of the fields, namely a = v and tt = 0, and average term 
by term. 

oo / \k —271 \ 

(F(a,7r))= E ^^'='"^^^,0) ^^) (17) 



Here 

jSjk+n 

We then need to relate {A"'tt^) to powers of (A^) and (tt^). Such relations 
were derived in [H]. We have (A") = for odd n and (A") = (r2-l)!!(A2)"/2 
for even n. For the pions, all the species are equivalent, so that (tt^) = (7r|) = 
{ttD = i(7r2) and (tt^'^) = (2A; + l)!!(|7r2)^ After substituting back into ([TTD, 
one notices that the averaging is equivalent to an integration over a Gaussian 
distribution, such that 

oo oo 

(F(a,7r2))= J dzP^iz) J dyy^P^{y)F{v + z,y^) (19) 

-oo 

where 

P^{z) = _— - exp ( -— — ) (20) 





PAy) = \- 7-^ -7T7±^ (21) 



Note that in the limit of vanishing mean square fluctuations {F{a,TT^)) 
F{v,0), as one would expect. 
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It is important to know that the approximation employed here is ther- 
modynamically self-consistent. There are two consistency relations between 
the mesonic masses and fluctuations that follow directly, namely 



oo 



{A^) = 2^ = ^ fdp^ J (22) 







oo 



and 



It is quite helpful to make use of a generalization of an expression in [12] for 
the derivative of f[T^ with respect to some parameter a. After two integra- 
tions by parts, one obtains 

We can use this to calculate the derivative of the thermodynamic potential 
with respect to the mean field as well as the sigma and pion quasi-particle 
masses. This gives 



2 



cff 



1 d{7i^) 1 




2 dv 






= 







2?l^i\^)-'"-H0 (25) 

These three derivatives vanish because of the averaged equation of motion 
and the two effective mass equations. This gives us the thermodynamic 
consistency we require. 

Due to thermodynamic consistency, we can calculate the energy density 
in the standard way, yielding 



{U) - \ml{A') - \ml{^') + dpp' 



6 7 , 2 / E E \ 

+-2jdpp {^^iEq;y^^ + ^mWTi/- ^^^^ 
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The baryon number density can be calculated as 





Note that in the quark terms there are pairs of angular brackets indicating 
the averaging procedure of eq. (ITTl) . This means that the momentum inte- 
grals are evaluated for a particular configuration of the mesonic fields with 
quark mass given by eq. ([7j), and then the mesonic fields are averaged over. 
One may choose to define a thermally averaged effective quark mass by the 
ensemble average of eq. ([7]), to wit 

m'^ml = / {{a') + (tt^)) = {v' + (A^) + (tt^)) . (28) 

However, the quark integrals are not evaluated with this single mass. 



4 Numerical Results 

The equation of motion and the effective meson mass equations must be 
solved self-consistently. The integrals in the thermal averages are calculated 
using an adaptive n-point rule algorithm, and the thermodynamic integrals 
are calculated using an adaptive Simpson's rule algorithm. We have verified 
numerically in a representative sample of combinations of T and fi that the 
thermodynamic identities 

6(/., T) = -P(/i, T) + T^^^^ + ^^nBi^^, T) (29) 

and 

nB{^^,T) = — (30) 

are satisfied. For more details on the mathematical derivations and numerical 
procedures see [2T] . 

To determine the location of a phase transition, one computes the pres- 
sure as a function of temperature and chemical potential. An example is 
shown in Fig. 1, where P is plotted as a function of fi for two different tem- 
peratures. In this example the vacuum pion mass is fixed at its physically 
observed value. For the higher temperature of 80 MeV there is only one 
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self-consistent solution to the equations described above. For the lower tem- 
perature of 50 MeV, there is a unique solution at large /i and another unique 
solution at small /x. For a range of fi centered about 900 MeV there are three 
solutions: one is associated with a continuation of the low density phase, one 
is associated with a continuation with the high density phase, while the third 
solution (not shown in the figure) is an unstable phase. The point where the 
two curves cross is the location of the phase transition. In this example it 
is first order since there is a discontinuity in the slope dP{^,T)/dfi = ub- 
Where each curve terminates is the limit of metastability for that phase. 
The thermodynamically favored phase is the one with the largest pressure. 
At some temperature between 50 and 80 MeV the slopes are equal at the 
crossing point, there are no metastable phases, and the second derivative 
is discontinuous. This is indicative of a second order phase transtion. The 
location of this point in the /i — T plane is the critical point. The critical 
point has coordinates designated {fic,Tc). 

In Fig. 2 is plotted the energy density divided by versus T, a com- 
monly used plot in studies of the equation of state at zero baryon chemical 
potential, for vacuum pion masses of and 200 MeV. For zero pion mass the 
phase transition is first order, whereas for a vacuum pion mass of 200 MeV 
there is no thermodynamic phase transition, only a crossover from one phase 
to the other. The contribution to the energy density for a gas of massless 
bosons is (7r^/30)T^ per degree of freedom, whereas massless fermions con- 
tribute 7/8 of that amount due to the difference between Bose and Fermi 
statistics. If all the degrees of freedom in this model were massless at high 
temperature the ratio e/T^ would be about 8.22. The results shown indicate 
a value of about 7.5 around a temperature of 200 MeV and not changing very 
quickly with increasing temperature. Obviously this is due to interactions, in 
particular the quasi-particles acquiring effective masses at high temperature 
due to thermal fluctuations. It should be remarked that if one considered a 
gas of massless quarks with 24 fermionic degrees of freedom (same as in this 
model) plus massless gluons with 16 bosonic degrees of freedom (versus 4 
in this model) one would have e/T^ approximately 12.17. Calculations with 
lattice QCD just above T = 200 MeV yield a value about 3/4 as large [22]. 
This indicates that the number of effective massless degrees of freedom is 
less, due to interactions, which is one argument for referring to the matter in 
this temperature range as strongly-coupled quark-gluon plasma. The sigma 
model with quarks exhibits a much smaller decrease in the effective number 
of degrees of freedom, suggesting that the fewer degrees of freedom incorpo- 
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rated in the sigma and pion fields versus the gluon fields may be qualitatively 
reasonable. 

We should also remark that when the vacuum pion mass is set to 138 
MeV (curve not shown) we reproduce the numerical results of reference |T6] , 
a check on our (and their) numerical work. 

The phase diagram in the /i — T plane for a sampling of vacuum pion 
masses is shown in Fig. 3. For the sake of discussion, pick one, for example 
'^TT.vac = 138 MeV. There is a curve of first order phase transition starting on 
the /i axis and arching to the left. This curve terminates at a critical point of 
approximately Tc = 75 MeV, /Xc = 850 MeV. Representative curves for m7r,vac 
= 50 and 200 MeV are also shown. For vac = there is no critical point; 
the curve smoothly extends from the /i to the T axis. For ^ac = 321 MeV 
the critical point sits on the T = axis, and for m7r,vac > 321 MeV there 
is no phase transition at all. Of course the precise numbers depend on the 
constants in this model, such as the vacuum sigma mass and what value one 
assigns to the coupling of the quark field to the sigma field, but the results are 
in line with expectations [12] . In particular, reference [3] found = 99 MeV 
and /ic = 621 MeV using the physically observed pion mass. Reference |3] 
used the mean field approximation whereas the present work also includes 
thermal fiuctuations in the meson and quark fields, and in addition that 
paper fixed the vacuum sigma mass at 600 MeV compared to the present 
value of 700 MeV. The increase in and decrease in fic compared to our 
results is consistent with a sort of universal curve, as exemplified by Fig. 3. 

Something very interesting happens as the vacuum pion mass is decreased 
from 50 MeV to 0. Figure 4 shows the phase diagram for a vacuum pion mass 
of 35 MeV. There are now two critical points! There is a line of first order 
phase transition beginning on the /x-axis and arching to the left to end at a 
critical point of fici ~ 725 Mev and T^i ~ 92 MeV, and another line of first 
order phase transition beginning on the T-axis and arching to the right to 
end at a critical point of fic2 ~ 240 Mev and Tc2 ~ 137 MeV. In the vicinity 
of /X = 575 MeV and T = 110 MeV the latent heat gets pinched to zero for 
some critical value of the vacuum pion mass between and 35 MeV. We have 
not attempted to pinpoint the exact value of the vacuum pion mass since it 
is quite numerically intensive and depends on all the other parameters in the 
model, as well as the structure of the model itself. 

The latent heat along the critical curve is plotted against T in Fig. 5 
for vacuum pion masses of 0, 35, and 138 MeV. At T = 0, corresponding to 
large chemical potential, the latent heats are on the order of 270 MeV/fm^. 
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Where they go to zero is a critical point. Note the nonzero latent heat in the 
bottom right corner of the figure; it corresponds to the region in the upper 
left corner of Fig. 4. In that region for that value of vacuum pion mass the 
latent heat is very small, less than about 10 MeV/fm^. For not much larger 
values of the vacuum pion mass this latent heat shrinks to zero. 

The results shown in Fig. 3 illustrate the conventional view on the nature 
and location of the critical point. The result shown in Fig. 4 is unconven- 
tional or exotic. To the best of our knowledge, this is the first time a model 
calculation has resulted in such a phase diagram. It is well known that the 
nature of the phase transition or crossover for two flavors of quarks is sensi- 
tive to such details as the strength of the axial U(l) anomaly and the value of 
the vacuum o mass. For example, see [221 EH [25] • It cannot be expected that 
a model as simple as the one studied here can make definitive predictions for 
what actually happens in QCD with the observed values of the quark masses. 
However, it can serve to illustrate the possibilities. 

5 Conclusion 

We have utilized a relatively simple yet self-consistent model to study the 
curve of chiral phase transitions in the plane of temperature and baryon 
chemical potential as a function of the vacuum pion mass. For a range of 
vacuum pion mass from zero to some critical value, which for the parameters 
we have chosen in this paper is less than 35 MeV, there is a continuous curve 
of first order phase transition from the T axis to the /i axis. As the vacuum 
pion mass increases, two critical points emerge; this is an unconventional, or 
exotic picture. Further increase in the pion mass causes one critical point to 
hit the T axis and disappear while the other critical point moves to smaller 
T and larger /i. For a pion mass of 321 MeV this remaining critical point 
hits the [I axis and also disappears. Changing the other parameters in the 
model, such as the vacuum a mass or the strength of the coupling of the 
quarks to the chiral fields, will obviously change the particular values of the 
pion masses for which two, one or no critical points exist. The existence and 
behavior of the critical point at large T and small is consistent with the 
lattice calculations of de Forcrand and Philipsen [lOj. (A recent preprint by 
Fukushima [26] shows that the Nambu Jona-Lasinio model with sufficiently 
strong vector coupling also supports a critical point at large T and small /z.) 
The existence and behavior of the critical point at small T and large [i is the 
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conventional picture. The model presented here shows that both pictures can 
coexist for a range of pion masses. More details on the numerical solutions 
as well as results when the aforementioned parameters are varied will be 
presented elsewhere [21] . 

Further work with this and related models are easily identified. It would 
be interesting to extend the sigma model to include strangeness, and to vary 
both the vacuum pion and kaon mass, using the self-consistent techniques of 
the present paper. It should be straightforward to calculate shear and bulk 
viscosities and thermal conductivity along the lines of [2^ and [28]. These 
transport coefficients could be used in fluid dynamic modeling of heavy ion 
collisions at FAIR and for low energy runs at RHIC. These projects are under 
current investigation. 
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Appendix 

To solve the equations of motion one needs the derivative of 

t^quark = — In Zquark)/!^ 

with respect to A 



quark 



9A 



g^{aA{m)) 



(31) 



(32) 



where the function A is the fermionic analogue of the mesonic mean square 
fluctuations 



A{m) = 2 



t^^quark 6 

dm? 



2 r 



+ 



(33) 



Within the current approach the quark condensate can be expressed as 

1 / <9?7quark 



^ \ da 

One also needs the second derivatives 



g{aA{m)) . 



(34) 



' FfiTT 
*^ '-^ quark 



g' ( A{m) + 2g^a^ 



OA- 



dA 



dm? 



2^ (^(" + A)/l(m)> 



' fi2TT \ 
^ ^ quark \ 



(A2 



(7rM(m)) (35) 



where the right hand sides are derived using the Hartree approximation in 
the scries representation of the thermal averaging process. The required 
derivative is 



dA 
dm^ 



3 °° 1 



E 



+ 



(36) 



One also requires the averages of the mesonic potential and its derivatives. 
For the pressure and energy density one needs the averaged potential itself. 



(37) 
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For the equation of motion one needs the first derivative of the bare potential 

(11)=^" + ^^^'^ + ^^'^ - / V ^ ' (38) 
and for the mesonic masses one needs the second derivatives 

(1^) = \{3v' + 3{A') + {^')-f) (39) 
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400 600 800 1000 

/2(MeV) 

Figure 1: Pressure versus baryon chemical potential for two different tem- 
peratures with the physical value of the pion mass in vacuum. At the lower 
temperature there is a discontinuity in the slope at the crossing point, indi- 
cating a first order phase transition. The curves are extended to the limits of 
met ast ability. At the higher temperature the pressure is smooth and there 
is no phase transition. 
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Figure 2: The energy density, normalized to T'^, versus T at zero chemical 
potential. For zero vacuum pion mass there is a first order transition, whereas 
for a vacuum pion mass of 200 MeV there is only a rapid crossover. 
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150 




/x(MeV) 

Figure 3: Phase diagram of this model for illustrative values of the vacuum 
pion mass. The solid curves represent a first order phase transition, ter- 
minating at a critical point where the transition is second order. For zero 
vacuum pion mass the curve is first order all the way from one axis to the 
other. 
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Figure 4: Phase diagram of this model for a vacuum pion mass of 35 MeV. 
The sohd curves represent a first order phase transition, terminating at a 
critical point where the transition is second order. For this value of the 
vacuum pion mass there are two critical points. 
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Figure 5: Latent heat along the critical curve for various vacuum pion masses. 
The critical point corresponds to the vanishing of the latent heat. Note that 
for a vacuum pion mass of 35 MeV there is a second region in the bottom 
right corner. 
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